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Foreword 


NASA-Langley’s Computational Materials program began in October 1994. The 
interdisciplinary team representing NASA Centers, universities, and industry met 
formally for the first time in January 1995 at LaRC. 

One year later, with some changes in personnel, the group assembled again, this time 
at The College of William and Mary, to review the technical progress that had been 
made and to identify future directions. The present document is a collection of the 
technical charts presented at that two-day workshop. 

The concept of Computational Materials requires model development at size scales 
ranging from atomic dimensions to the macroscopic, and times ranging from 
picoseconds to weeks. The LaRC program therefore adopts a hierarchical approach, 
with careful attention paid to the links between disciplines. This collection of papers is 
presented in order of increasing scale, as follows: 



Discipline 

Scale 

G. D. Smith, 

U. Missouri-Columbia 

Quantum/molecular 

mechanics 

Subatomic to molecular 

B. L. Farmer 
U. Virginia 

Molecular mechanics 

Molecular 

E. J. Siochi 
Lockheed 

Physical chemistry 

Molecular 

J. A. Young 
U. Virginia 

Molecular dynamics 

Molecular 

assemblies 

E. J. Dawnkaski 
College of William & Mary 

Monte Carlo 
simulation 

Molecular 

assemblies 

R. A. Orwoll 

College of William & Mary 

Chemical Thermodynamics 

Liquids 

T. Yost, J. Cantrell 
NASA LaRC 

Physics 

Solids 

W. W. Welsh 

Chemistry; statistics 

Molecular to continuum 

W. G. Knauss 
California Inst, of Tech. 

Constitutive 

models 

Continuum 

T. S. Gates 
NASA LaRC 

Viscoplasticity 

Micromechanics to 
macroscopic 
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Structure of the model 
ether-imide 
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Bonded parameters 
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Bonded parameters, c 
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T or sional char acteri Stic s 
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Partial atomic charges 
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Partial atomic charges, 
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Nonbonded parameters 
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Nonbonded parameters, cont. 
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Nonbonded parameters 
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Out-of-plane deformation 
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Molecular dynamics simulations 
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Density 
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Thermal expansion 
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Sensitivity of density to force 

field parameters 
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Pair distribution functions 



60 

C 

• fH 

a 

e 

• 

03 

ts 

a 

ctf 

o 

X 

& 

<D 

a 

>> 


X) 


n ^ 

o S 

too .y O 

g o 

2 g g 

o g § 

-Sc 15 

— , -H 

5 H O 

S)^ ~ 

<d 3 g 

o O i 

V ^ 

> O 05 

JS g ^ 

<u S .g 
^18 


19 


University of Missouri-Columbia 4 January 


Pair distribution functions, cont 
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Pair distribution functions, cont. 
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Ongoing work 
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Key issues and future work 
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Table I 


Quantum Chemistry Based Force Field for Ether-Imide 


Type 

Energy Function 

Parameters 

Bond Stretch 

Constrained 


(A) 


C—O (carbonyl) 3 



1.19 


C-0 (ether) 



1.38 


C-H 



1.07 


C=C 



1.38 


CfN 



1.43 


c-c, 



1.49 


Valence Bend 

E ijk = 0.5 ky k (0 yi - e o yi ) 2 

k, Jk (kcal/mol rad 2 ) 

o°«. 

c-c r N 


144. 


1.817 

c-c-c, b 


100.8 


1.882 

c-c-q c 


100.8 


2.286 

Cj-N-Cj 


115.2 


1.988 

c-c.-o 


144. 


2.253 

c-c-c 


144. 


2.094 

C-C-H 


72. 


2.094 

C.-N-C 


115.2 


2.148 

C-O-C 


148.6 


2.072 

C-C-0 


100.8 


2.094 

N-Cj-O 


100.8 


2.1991 

N-C-C 


115.2 


2.094 


25 






Torsion E ijkl = 0.5[k 2 ijkl (l-cos(2<J> ijU )) + k 4 ijtl (l-cos(4<t» ijn ))J 


k 2 1Jtl (kcal/mol) k 4 1Jld 


c-c-c-c 

25.92 

0.00 

C-C-C-H 

0.00 

0.00 

H-C-C-H 

0.00 

0.00 

C r C-C-C 

0.00 

0.00 

C-C-C-H 

0.00 

0.00 

Q-C-C-C, 

0.00 

0.00 

o-c r c-c 

0.00 

0.00 

N-C-C-C 

0.72 

0.00 

c-N-c r c 

0.00 

0.00 

C-N-Cj-0 

14.40 

0.00 

C-C-O-C 

0.86 

0.00 

C-C-C-0 

0.00 

0.00 

H-C-C-O 

0.00 

0.00 

C-C.-C-O 

0.00 

0.00 

C r N-C r O 

7.06 

0.00 

c-c-N-q 

1.20 

0.26 

N-C-C-C 

0.00 

0.00 

N-C-C-H 

0.00 

0.00 
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William & Mary Computational Materials Workshop 
January, 1996 
Barry L. Farmer 
Retrospective Summary 

(including both material presented and insight gained from resulting discussions) 

The UVA (on-grounds) effort has focused on exploring the validity of the Tripos force field for 
computing the interactions between a graphite surface and a polyimide chain. Validation has proceeded along 
two paths. 

The first validation path has involved comparison of calculations using the Tripos force field with the 
results of semiempirical molecular orbital and ab initio calculations for pairs of molecules, e.g. methane, 
formaldehyde, and trimethyl amine with benzene. The comparisons have not been good. However, rather than 
point to a deficiency in the empirical (Tripos) force field, the results underscore the inability of the molecular 
orbital and ab initio methods to properly represent dispersion interactions except when very extensive basis sets 
are employed. Thus the validation of the Tripos force field lies more correctly in the comparisons with 
experimental data than with the results with the more fundamental computational methods. 

The second validation pathway involves calculations of the association between various small molecules 
(corresponding to representative fragments of polyimides) and a section of a graphitic surface (specifically, a 
coronene molecule) were made. These were compared with previous computational results based on more 
elegant treatments of the interaction potential developed to reproduce experimental data. Both the interaction 
energy versus distance curve and the geometry of the minimum energy association calculated using the Tripos 
force field for benzene on graphite were in good agreement with the previous model and data. In addition, 
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comparisons of the experimental heats of adsorption of other molecular fragments with graphitic substrates were 
generally in good agreement with the calculated interaction energies. 

Another aspect of the semiempirical molecular orbital and ab initio calculations results however is both 
useful and interesting. Graphite, as a material, is electrically conductive. The effects of this conductivity on 
graphite’s interactions with adsorbing material has apparently not been studied previously in the context of 
modem computational methods. (This neglect may arise in fact from the inability of molecular orbital 
calculations to provide a reasonable treatment of dispersion forces). A dipole in the vicinity of a conducting 
(metal) surface interacts electrostatically with an image dipole as well as through dispersion effects with the 
substrate atoms. The image interaction is about an order of magnitude smaller than the dispersive ones. This 
indicates that, overall, “summation of dispersion and repulsive energies” may indeed be a good approximation for 
computing the interaction energies of adsorbate molecules with graphite. At the same time, ignoring the image 
interactions may lead to underestimation of the role of the dipolar moieties of the polyimides on the interactions 
with the graphite substrate. The size of these effects will determine whether we will need to modify the potential 
function parameters to get a reasonable representation of the behavior of a polyimide on a graphite substrate 
Since the interaction is electrostatic, the results of semiempirical molecular orbital and ab initio calculations 
should be valid in spite of the inability of these methods to properly account for dispersion interactions. 

The calculations to date indicate that the partial charges in the adsorbing molecule and those on the atoms 
of the “substrate” molecule (benzene, naphthalene, coronene) change to a small extent as the separation distance 
changes. As expected, the charges on the substrate atoms nearest the adsorbate change individually more than 
those at greater distances, drawing upon contributions from its own several neighbors. The perturbations in 
electron distribution extend out from the location of the adsorbate, but they are damped in magnitude. 
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Figure Captions and Annotations 


Note: 

Sybyl does molecular mechanics calculations using a classical description of the force field. 

Spartan does semiempirical molecular orbital (quantum mechanical) calculations. 

1. Title 

2. Polyimides of Interest 

3. Energy (and partial atomic charge on methane hydrogen) versus separation distance between benzene and 

methane calculated using Spartan (semiempirical molecular orbital methods). 

4. Energy (and partial atomic charge on fluorine) versus separation distance between benzene and 

perfluoromethane calculated using Spartan (semiempirical molecular orbital methods). 

5. Energy (and partial atomic charge on nitrogen) versus separation distance between benzene and 

thrimethylamine calculated using Spartan (semiempirical molecular orbital methods). 

6. Energy (and partial atomic charge on oxygen) versus separation distance between benzene and formaldehyde 

calculated using Spartan (semiempirical molecular orbital methods). 

7. Charges on the formaldehyde carbon and oxygen and three of the benzene carbons as a function of distance 

between benzene and formaldehyde (calculated using Spartan semiempirical molecular orbital methods). 

These date show that the proximity of a polar (carbonyl) group perturbs the distribution of charges in the 
benzene molecule. What might be the effect of polar groups proximate to delocalized ring systems such as 
naphthalene, anthracene, coronene, and ultimately, graphite. 

8. Charges on the formaldehyde carbon and oxygen and three of the naphthalene carbons in Ring 1 as a function 

of distance between naphthalene and formaldehyde (calculated using Spartan semiempirical molecular 
orbital methods). The formaldehyde molecule was centered over Ring 1 of the naphthalene molecule. 

9. Charges on the formaldehyde carbon and oxygen and four of the naphthalene carbons in Ring 2 as a function 

of distance between naphthalene and formaldehyde (calculated using Spartan semiempirical molecular 
orbital methods). The formaldehyde molecule was centered over Ring 1 of the naphthalene molecule. 

10. Charges on the formaldehyde carbon and oxygen and three of the anthracene carbons in Ring 1 as a function 

of distance between anthracene and formaldehyde (calculated using Spartan semiempirical molecular 
orbital methods). The formaldehyde molecule was centered over Ring 1 of the anthracene molecule. 

11. Charges on the formaldehyde carbon and oxygen and three of the anthracene carbons in Ring 1 as a function 

of distance between anthracene and formaldehyde (calculated using Spartan semiempirical molecular 
orbital methods). The formaldehyde molecule was centered over Ring 2 of the anthracene molecule. 
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12. Charges on the formaldehyde carbon and oxygen and four of the anthracene carbons in Ring 2 as a function 

of distance between anthracene and formaldehyde (calculated using Spartan semiempirical molecular 
orbital methods). The formaldehyde molecule was centered over Ring 3 of the anthracene molecule. 

These data show that a polar group does induce charge redistribution in rings sharing delocalized electronic 
structure with the ring closest to the approaching polar functionality. The effect, not surprisingly, dies out with 
distance. This can be attributed (in part if not totally) to the fact that there is an increasing number of atoms to 
contribute decreasing amounts of charge (per atom). 

How do molecular mechanics (Sybyl) results compare with semiempirical (Spartan) results. Specifically, how 
do the energy versus distance curves compare? 

13. Energy versus separation distance between coronene and methane calculated using Tripos force field and 

Sybyl for molecular mechanics. 

14. Energy versus separation distance between coronene and methane calculated using Spartan semiempirical 

molecular orbital methods. 


15. Energy versus separation distance between coronene and formaldehyde calculated using Tripos force field 

and Sybyl for molecular mechanics. 

16. Energy versus separation distance between coronene and formaldehyde calculated using Spartan 

semiempirical molecular orbital methods. 

The molecular mechanics results and semiempirical molecular orbital results are not in good agreement. 
Molecular mechanics methods compute van der Wools non-bonded interactions and Coulombic interactions as 
a matter of course. Semiempirical molecular orbital methods can account mainly for the Coulombic effects. 
Similarly, even ab initio methods cannot take into account non-bonded interactions (because they arise from 
electron correlation effects) unless prohibitively large, delocalized basis sets are used. 

Comparisons with experimental data and other types of calculations in the literature might provide a better 
measure of the validity of the Tripos force field 

17. Orthogonal views of a benzene molecule proximate to a coronene molecule (representing a section of a 

graphite sheet). 

18. Relative energy as a function of distance between benzene and coronene (in the positions shown in the 

previous figure). 

19. Title and figures from a paper reporting calculations (and comparisons with experimental data) for the 

absorption of benzene on graphite. 

The Sybyl results are in very nice agreement with the results in this paper, both in regard to the minimum 
energy orientation, the depth of the energy minimum, and the minimum energy separation distance. 
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Molecular Mechanics 

and 

Molecular Dynamics 

with applications to polyimides as 
piezoelectric and composite materials 

B. L. Farmer 

Materials Science and Engineering 
University of Virginia 
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Piezoelectric Polyimides 



32 




33 


Distance (angstroms) 


Energy and Charge Versus Distance 
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Energy and Charge Versus Distance 

Phen_nndmn Molecule 
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OUTLINE 
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Polyimides as materials of choice 
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Universal Calibration Curve 
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Effect of Molecular Weight on the 
Properties of a Soluble Polyimide 
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Find relationships between physical properties of polymers and 
molecular weights 



Differential Scanning Calorimetry 
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Profile of Vacuum Microextruder 



64 


Single Screw 



Melt Viscosity Calculation 
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= Capillary Length 



Relationship of Melt Viscosity to Molecular Weight 
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Microhardness Determination 
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p = Load (kg) 

AP = Area of the Impression 



Relationship of Hardness to Molecular Weight 
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Density Determination Apparatus 
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= Weight of Solid Body in Air 
= Buoyancy of Solid Body in Test Liquid 
= Density of Test Liquid at a Given Temperature 
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Compact Tension Configuration for Fracture Toughness 
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Effect of Stoichiometric Imbalance on 
Molecular Weight Distribution 
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SEM Photomicrographs of LaRC™-SI Fracture Surfaces 
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Fracture toughness decrease) 
offset in stoichiometry. 
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Molecular Modeling of Novel Piezoelectric 
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PIEZOELECTRICITY 
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* Joycelyn 0. Simpson et al., Novel Piezoelectric Polvimides . Proceedings of the 1995 Fall Meeting 
of the Materials Research Society (Novemer27-December 1, 1995) Boston, Massachusetts. 







PARAMETERIZA TION 


co 

G 

o 

* 

S 

3 

_o 

3 

p 


2 

13 


P 

P 

£ 

S3 

8 

a. 


p 

“ 4-» 

0) ^ 

SS co 
cd o 

■° c 

2 I 

2 b 

P ^ 

8 c3 

2 3 

w p 

s 2 

o o t 

b £ C 

§ E 5 

p 

o o o 


o s 

CO X 


a> 


C P 

_ S x> 

X X +-* 

cd O co 

§? £ £ 

*53 b 

x> «* 


° 3 


<D 


v. 2 

O o ^ 

A c £ a 

Dh w a 


I 

cd 

u 


CO 

g o 

X _c2 

w t3 
4-h P 

° 6 

CZ) . 

c g 

a .8 

G c 

O S • 

G-( «-G —j 

G p 2 

o 8 ® 

O B 

*S I 
m 3 o 
8 

J S ” 

S .2 -rH 

Pf 'G _C2 
« 2 .i a 
G cd £■ 
P 


P 

P 


cd 


cd 


G 

o 


CO 


_ cd *G 
G P § 
2 P 43 

0 G P 
Cl O P 

~7 X G 

0> k> c 

*1 

P X 

J— i 

ed P 

01 T3 

E c 

o .5 

U >» 


3S 


Vh 

2 
3 

p 

2 

3 

E 

5* T3 
O t3 

Oh Cd 


2 

3 


p 

P CO 

b p 

tQ u 

*** ,cd 

ON p 

□Li co 

E o 

U x 
« 

CO f? 


00 

O 

HH 

PQ 

«4i 

o 

G 

O 

• 

W 

cd 

o 

• f-H 
•w* 

T3 

O 


❖ 


X) 

< 

U 

< 

Gh 

o 


p 


*0 

P 

P 

G 

T3 

O 

Vh 

P 

t-i 


82 





CD 


a, cs 



83 







03 

cd 

£ 

X 

o 

x 

c« 

3 

o 

& 

§g 

Cd 

o CO 

5g O 

•c ® 

^ c 

Oh •’""' 
<D 


© 

73 

o 

s 

13 

u 


cd 

c 

o 

• fH 

C/3 

a 

E 

'-5 g 
o o 
5 x 

f & 

cd o 
<*5 S 

c < 
.1 £ 


o 


cd 


2'? 

3 g 

® § 

U 73 


<D 
O 

0) 

T3 

• fH 

s 

•fH 

>* „ 

T3 © 


O 
Oh 

Wh 

8 

a> 

> 

<4H 

<rP 

s 

o 

"£b 

CN 

O 


03 

g 

CD 

73 

£ 

0 

1 < 

cd 

2 

o 

PQ 


9£ ae 

ww uu 


c 

o 

C/3 

o • 

E^p 
cd G 
G ^ 

^ W> 

VH ^ 

cd cn 

3 w 

0 w 
u >» 

g-3 

1 § 

X cd 

O G 
O <D 
O G 

m .5 

W" u- 

<D £ 
J-h Oh 

G X 

03 <D 

C/3 

CD 03 

js 2 

.» 3 

X *-« 

§ B 

^ CD 


03 

>> 

03 


°.g 

C/5 +-> 

« <D 

.cd t8 


■ E 

X .53 

O “ 

■§1 

C C/5 

2 s 

cd 


>> Qj 

5 « 

G C/3 
G a) 

2 x 

73 § 
M co 

O D 

cd &0 

^ S3 

*G 

73 O 

-? n W) 
Ot) c 
\G -S 

2 a 

o 2 

03 
CD 

cd 33 

U 8 

W) CD 

6 E 

2 ^ 

o 

Oh 


3S 


84 


electric poling field. 



Platin 
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View of parent box (red) within the systems periodic boundaries. 
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fhe parent cell containing a five-mer polyimide and plating atoms. 


96Molecular dynamics was done to obtain a starting configuration. 
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POLE MATERIAL 
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Polarization of Cell 
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A small polarization of 0.006 C/m2 was obtained when using 

experimental parameters. 


Collapse Time Scale 
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Dipole Reorientation 
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Z Component 



Polarization of Cell 



Time (ps) 

A polarization of 0.02C/m2 was obtained when using the collapsed 

time scale oarameters 



CONCLUSIONS 


o 


"O ^ 

s g 

O *1 
OJ — 

F cd 

£ o 

0 

° § 

a .2 

v o 

* 1 

<Z> ** 

S I 

1 I 
§ §■ 


^ X 

X) J=3 

II p 

TD o 

X CO 

P p 
T 3 O 
o X 

s e- 

CO O 

cd C 

* I 

<D o 

3 '-3 
S .2 

•r* *-< 
>* P 
X Oh 


C 

P 

<D 

X) «5 
P V 

> U 

c« *a 

X G 

• f*H 

CO 

P X< 
S 3 O 
O" Oh 

G 22 

X 5 

o o 

<D X 


cd Cd 
X X 

s ° 

X bo 

co G 


5P « 

G <d 

S s 


*G cd 

S -S 

cd (D 

CO ^ 

P bo 

b a 
3 o 

cd ^ 

S _g 

O-i X 


T3 X 

a. ^ 

^ X 

P 

P £ 

;*“H 


P o 

p Oh 

H— > 

X 

<-t-H H— * 

p § 

a I 

t£ & 

G g 

^ .2 

0) oo 

H | 

Oh P 

6 -S 

o _ 
O « 

w> 
• *— 1 

55 

X 

0 & 

©1 

0 1 



95 


K Polarizations of 0.02 C/m z were calculated for the 
polyimide. 



FUTURE WORK 





96 


piezoelectric responses, calculation ot absolute values is no 
focus of this research. 
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Monte Carlo Similation of a Polyimide Melt 


http : //pooh. chem. wm. edu/~ej d/resear ch/meetl/ 
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01 / 05/96 10 : 49:47 



Why Monte Carlo? 


http://pooh.chem.wm.edu/. . . esearch/meetl/pagel .html 


Monte Carlo 

■ Bridge the gap between sub-microscopic and macroscopic 

□ Retain knowledge and detail of sub-microscopic system 

□ Determine macroscopic quantities 

■ Larger system sizes 

■ Faster simulations (larger time frame available) 
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01/05/96 10:49:59 



What Type of Monte Carlo? (A) 


http://pooh.chem.wm.edu/. . . esearch/meetl/page2 .html 


Types of Monte Carlo 


■ Slithering Snake 

■ Kink-Jump 

■ Bond Fluctuation 
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Slithering Snake Monte Carlo 


http: //pooh. chem.wm. edu/ . . . /meetl/slithersnake .html 


Slithering Snake Monte Carlo 



■ One move type 

■ Unrealistic on local scale 

T. Wall and F. Mandel, J. Chem. Phys. 63, 4592 (1975) 
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Kink Jump Monte Carlo 


http://pooh.chem.wm.edu/. . . arch/meet II kink jump .html 


Kink Jump Monte Carlo 


j/' kink-jump 


end-bond 


IZJ 




s t 


crank shaft 


H. Verdier and W. H. Stockmayer, J. Chem. Phys. 36, 227 (1962). 
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Bond Fluctuation Monte Carlo 


http://pooh.chem.wm.edu/. . . etl/bondf luctuation.html 


Bond Fluctuation Monte Carlo 



P. Deutsch and K. Binder, J. Chem. Phys. 94, 2294 (1991). 

Paul, K. Binder, J. Batoulis, B. Pittel and K. H. Sommer, Makromoi Chem., Macromol Symp. 65, 1 (1993). 
Schulz and J. U. Sommer, J. Chem. Phys. 96, 7102 (1992). 

4. J. U. Sommer, M. Schulz and H. L. Trautenberg, J. Chem. Phys. 98, 7515 (1993). 
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Bond Fluctuation MC vs Kink-Jump MC 


http://pooh.chem.wm.edu/. . . es ear ch/meetl /page5 .html 



BFMC 

KJMC 

ADVANTAGES | 

□ Computationally more efficient 

□ One move type 

□ Acceptance probability higher 

□ Faster equilibration in dense systems 

□ Longer move length 

□ Faster equilibration in less dense systems 

DISADVANTAGES | 

□ Short move length 

□ More steps to reach equilibrium 

□ Slower equilibration time in less dense 
systems 

□ Three types of motion 

□ Lower acceptance probabilities 

j 

□ Computationally less efficient 

□ Slower equilibration time in dense 
systems 
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More Bond Fluctuation Monte Carlo 


http://pooh.chem.wm.edu/. . . esearch/meetl/page9 .html 


Bond Fluctuation Monte Carlo: 

The details 

1 

! 

Allowed motion: | 

1 

! 

This yields 108 Bonds with 5 
different bond lengths 

P( 1,0,0) 1 

b = 4 1/2 (2), 5 m , 6 1/2 , 9 1/2 (3), 
101/2 

Allowed bond vectors: 
Union of: 

P(2,0,0), p(2,l,0), 
| p(2,l,l), p(2,2,l), 
j p(3,0,0),p(3,l,0) 

Two successive bonds may have 
up to 87 different angles over the 
interval [0-180] degrees. 
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Allowed Bond Angle Distribution 


http://pooh.chem.wm.edu/. . . search/meetl/pagelO . html 


Allowed Bond Angle Distribution 
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BFMC restrictions summary 


http://pooh.chem.wm.edu/. . . search/meetl/pagell .html 


Summary of Restrictions 

Excluded Volume: 

The surrounding 26 grid locations must be vacant. 

Bond Vector: 

Set of allowed bond vectors prevents intrachain crossing. 
Allows only 5 bond lengths. 

Allows only 87 bond angles. 

What About Energetics? 

Metropolis Transition Probability: w = min[l, 


107 


01 / 05/96 10 : 54:51 



BFMC energy restrictions summary 


http://pooh.chem.wm.edu/. . . search/meetl/page!3 .html 


Energy Restrictions 

□ Psuedo Bond Length 

No Restrictions 

□ Bond Angles 

Restricted by bond bend stiffness (calculated via biosym) 

□ "Long Range Torsions” 

Calculate Energies for torsions and restrict accordingly 

□ Inter-chain interactions: 

Not accounted for as of yet 


108 


01 / 05/96 10 : 55:10 



What is LARC-IA? 


http: / /pooh, chem.wm. edu/ . . . esearch/meetl/page6 .html 


Go to Page Seven 

LARC-IA 
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Energy (Kcal) 




Bond Fluctuation Representation of LARC-IA 


http: / /pooh, chem.wni. edu/ . . . esearch/meetl/page8 .html 


Bond Fluctuation representation of LARC-IA 
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What is LARC-IA? 


http://pooh.chem.wm.edu/. . ,esearch/meetl/page7 .html 


Six sub-units of one monomer of LARC-IA 
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Cross-linking 


http://pooh.chem.wm.edu/. . . search/meetl/pagel9 . html 


Cross Linking 


□ Allow only end cross linking. 

□ Chain Ends will be allowed to cross link with other chain ends. 

□ All restrictions must still be met. 

□ Cross link kinetics as function of MC step (and also Time?). 
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Computational Materials Workshop, January 4-5, 1996 


Equation of State Measurements: 
Thermal Pressure Coefficient 
of an Imide-Ketone Model Compound 

Robert Orwoll and Rachel Ward 
Department of Chemistry 
College of William and Mary 
Williamsburg, VA 23187 
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Thermal Pressure Coefficient = y v = — 
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THERMODYNAMIC RELATIONSHIPS CONNECTING THE DESIRED (8E/aV) T AND THE 
MEASURED THERMAL PRESSURE COEFFICIENT (8p/8T)v 
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Transistor Relay 



Tungsten Wire 


“CONSTANT” VOLUME CELL FOR MEASURING THE THERMAL PRESSURE 
COEFFICIENT 
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o 


MOLECULAR STRUCTURES FOR THE MODEL ETHER-IMIDE (TOP) AND KETONE- 
IMLDE (BOTTOM) 
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Pressure (psi) 


Pressure vs Temperature 

of a Model Imide-Ketone 



RESULTS FOR THE KETONE-IMEDE (BEFORE CORRECTION FOR THE EXPANSION 
AND COMPRESSION OF THE MERCURY AND GLASS) 
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Comparison with Other Amorphous Compounds 



(5E/dV) x , J/cm 3 

Temp, °C 

/i-octane 

265 

25 

benzene 

379 

25 

cyclohexanone 

413 

25 

poly(dimethyl siloxane) 

145 

180 

polyethylene 

275 

180 

polyethylene oxide 

433 

100 

imide ketone 

440 

186 


p-V-T Results to Date 


Model Cmpd 

CTE = (l/V)(dV/dT)„, K 1 

(aE/av^, j/cm 3 

imide ketone 

6.59X10 -4 (±3x10-®) 

440±10 

imide ether 

6.83X10" 4 (±3x1 O'®) 

next 

measurement 
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Relationship between 
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is the particle location (at rest) 

.is the stress acting on the segment 





Hooke's Law 
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co is the angular frequency of the wave 
(= 27cf , f is frequency) 

i is the square root of -1 



Elastic Constants and 
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Thermodynamic Definitions 
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Thermodynamic Definitions 
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Thermodynamic Definition 
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Viscoelastic Materials 
and Ultrasonic Measurements 
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velocity can be determined. 

This ultrasonic technique measures the unrelaxed second-order moduli of 
the LaRC-SI sample 
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Ultrasonic Shear Wave Velocity 

in LaRC-SI 
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Molecular Weight Offset 
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Isentropic Shear Wave Modulus 

in LaRC-SI 
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Molecular Weight Offset 




Ultrasonic Wave Compressional Velocity 

LaRC-SI 
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Molecular Weight Offset 




Isentropic Compressional Wave Modulus 

in LaRC-SI 
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Molecular Weight Offset 




Measured Density 
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Wave Velocities as a Function 
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Moduli from Wave Velocities 



ISO 


1 . 6341 ±.006 1 . 6358±.006 8 . 978+.03 




Moduli Derived from Ultrasonics 
Isentropic Bulk Modulus 
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6 . 798+.05 


Conclusions 


a> 


CO 


H- C 

or® 

Us 

c s CO 

CD 

■o^a 

1|S 

3 -o J- 

O Q> $3 
JZ N m- 
WZ O 

*D J2 <D 

COE 

<o a. g 

W Jr CO 
C®u 

.2? o 

» 3 O 

o 

.£ co co 

§s§ 

M ^ Qi 

2 co co 
^ — a> 
co « *- 
a> 3 co 

■5.2 •" 

E co .2 
o> sz 

W h. 4 -* 

O) 


a> ® -5 

^ co 2 

o o 

JS co E 

-5 o) 

Pf 

CO g o 

(0 

CTJC 

cd zz o> 

E co “ 

a) > a 

ir (0 O 
3 CD 
co .co & 

CO 

<D >» w 

E of 

co q 2 

<D § CO 
£« ® 

CO «*- i- 

b 

® J- 

C Q.~ CO 

■ E * S 

Q D C O 
° 0) <D O 
CD CO > v. 
■C CO CD Q. 


153 


moduli from ultrasonic measurements. 
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GOALS & OBJECTIVES 


DEVELOP AND APPLY METHODS IN COMPUTATIONAL 
CHEMISTRY FOR FAST, RELIABLE, PREDICTION OF 
POLYMER PROPERTIES 

RELATE BULK MACROSCOPIC POLYMER PROPERTIES TO 
MICROSCOPIC MOLECULAR STRUCTURE AND PROPERTIES 

DEVELOP STRUCTURE-PROPERTY RELATIONSHIPS 



COMPUTATIONAL APPROACHES 


QUANTITATIVE STRUCTURE-PROPERTY RELATIONSHIPS (QSPR) 

-4 Group-Additivity Techniques (van Krevelen, Dow) 

-4 Connectivity Indices (Bicerano’s Synthia) 


MULTIVARIATE REGRESSION 

-4 Partial Least Squares (PLS) 

ARTIFICIAL NEURAL NETWORKS (ANNs) 


ATOMISTIC MOLECULAR MODELING 


-4 Cerius 2 Mechanical Properties Module 



DESCRIPTION OF METHODS 
QSPR Approaches 
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QSPR Approaches 


The van Kreveien Method 

D. W. van Kreveien, Properties of Polymers: Their Estimation and Correlation with 

Chemical Structure (Elsevier, 1976) 

• assumes that the physiochemical properties of a polymer can 
be obtained by adding contributions from the constituent 
chemical groups 

• this approach is commonly known as Group Additivity (GA) 

• based on empirical and semi-empirical fitting to expt’l data 

• relies on large database containing group contributions, 
typically one for each property of interest (e.g., T g , a Tl ) 

• property values apply to bulk amorphous and semi-crystalline 
homopoioymers, statistical polymers, and for polymer 
solutions 

Advantages 

fast 

easy to use 

-> applicable to many properties 
usually reliable 

Disadvantages 

depends on chemical-fragment library 

_> theoretically less rigorous than atomistic methods 
do not account for special interactions (e.g., solvent) 
do not provide insight into mechanism 
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Table 9-12. QSPR Groups (Page 1 of 3) 
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Table 9-12. QSPR Groups (Page 2 of 3) 


N.</ 

N 'c-e'" 
h / x ch 3 

\T 

1 

A 

9 

u> 

N c-c / 
h / x ci 

IS.transethene 
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27 .transmetby letfaene 
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a 

V 

1 

U 
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-C- 

u 

0 

5 

29,traiisdiioroethene 

30.ethyne 

31. carbonyl 

32.onbopbe©yiene 



/C H3 

p 

Xr 


'CI% 

ol 
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35 ,26dimetbppbeny»ene 
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C3 3 
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■>»■ 
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Table 9-12. QSPR Groups (Page 3 of 3) 


N N 

it ii 

N N 

HC S 
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I! 

ii ii 

-c c- 

\ / 
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H H 
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\ / 
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II U 
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<a>- 

-<a> 

H 
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H H 
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56. benzodioxazole 

<K> 

ok> 

H 

ch 3 

-O-'O 

CH 3 


57. parabiphaxylene 

58. pbisphcnylmetfaylene 

59. pbispben yiA 
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OO 
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0^0 
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ch 3 
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6 
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<>K> 
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CHi CH 5^3 
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CH 3 CH CH 3 
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-CH- 
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ch 2 

CH 3 0*3 
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67. 27naphthaicae 

58. fbutyimerbykse 
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van Krevelen Methodology 

The polymer properties that can be calculated with the van Krevelen methodology arc: 

Thermophysical Properties 

• glass transition temperature 

• crystalline melt transition temperature 

• attainable degree of crys tallini ty 

• maximum linear growth rate of spent litic crystallites 

• density of glassy, rubbery and crystalline phases 

• volume coefficients of thermal expansion 

• volume change on melting 

• molar heat capacity at constant pressure 

• molar heat capacity at constant volume 

• molar entropy of fusion 

• molar latent heat of fusion 

• cohesive energy 

• solubility parameter 

• surface tension 

Mechanical Properties 

• bulk modulus 

• Poisson’s ratio 

• longitudinal velocity of sound 

• shear modulus of an amorphous or semicrystalline polymer 

• tensile modulus of an amorphous or semicrystalline polymer 
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Ultimate Mechanical Properties 

• tensile yield strength 

• compressive strength 

• flexural strength 

• indentation hardness 

Transport Properties 

Permeation of a Gas 

• specific permachor 

• permeability of a gas 

Viscoelastic Properties of a Polymer Melt 

• activation energy of viscous flow 

• Newtonian viscosity 

• characteristic deformation time 

• non-Newtonian viscosity 

Dilute Solution Properties in a theta Solvent 

• Maric-Houwink prefactor 

• intrinsic viscosity 

• radius of gyration 

• critical molecular for entanglement 

Dilute Solution Properties in a Good Solvent 

• Maric-Houwink prefactor and exponent 

• excluded volume expansion factor 

• intrinsic viscosity 

• radius of gyration 

Concentrated Solution Property in a Good Solvent 

• Newtonian viscosity 

Electrical, Optical, and Magnetic Properties 

• dielectric constant 

• resistivity 

• refractive index 

• magnetic susceptibility 
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Properties Relating to Thermal Stability 

• free enthalpy of formation 

temperature ai half thermal decomposition 

• char residue 

• oxygen index 


The van Krevelen methodology relies on the calculation of various molar properties of 

O The l g ® repeatUmL ^ejnolar properties are denoted here by enciosatSTSLs 
{>. The molar properties are then used in the calculation of the macroscopic properties 
o e polymer. The molar properties that are calculated are l isted below: 

• mass {M} 

• number of QSPR groups { N } 


• number of backbone atoms {2> 

• van der Waals volume {V w } 

• glass transition temperature {Y g } 

• melt transition temperature { Y m > 

volume at 298 K in the amorphous pha<a» {y a y 

• heat capacity of the solid {C prJ } 

• heat capacity of the liquid {C p j} 

• entropy of melting (AS m > 

• cohesive energy {E^} 


• interaction <F> 

• parachor {P s } 

• permachor {II > 


• elastic wave velocity { U } 

• intrinsic viscosity {J> 


viscosity-temperature gradient {H h } 
polarization {P^} 
optical refraction {R GD } 
magnetic susceptibility {%} 
enthalpy of formation {A C a } 
entropy of formation {AC b } 


half thermal decomposition temperature {Y^ yj} 
char forming tendency {C^} 
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QSPR Approaches 


The Dow Method 

J. Seitz, J. Appl. Poiym. Sci., 49, 1331 (1993) 


• adopts group-additivity (GA) concept, but uses only six group 
contributions of the repeat unit to predict polymer properties: 
(1) molecular weight; (2) length; (3) vdW volume; (4) cohesive 
energy; (5) rotational df of the backbone; and (6) T g 

• focuses on mechanical properties 

• limited to amorphous (not semi-crystalline) polymers 


Advantages 

-» same advantages as van Krevelen: fast, reliable 
-> less dependent on large database of group contributions 
-» better fitted to high-performance engineering polymers 
-» conceptually more rigorous than van Krevelen 

Disadvantages 

-» same disadvantages as van Krevelen 
restricted to amorphous polymers 
-> depends on library of group contributions, although smaller 
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Properties Calculable 


Dow Methodology 

^5 ^"methodology implemented within the QSPR module it is possible to 
pre c different properties of amorphous thermoplastic polymers. These are: 

• Molar volume per repeat unit 

• Density 

• Thermal expansion coefficient 

• Cohesive energy 

• Solubility parameter 

• Surface tension 

• Dielectric constant 

• Average molecular cross-sectional area 

• Glass transition temperature 

• Poisson’s ratio 

• Tensile modulus 


• Tensile yield strength 

• Brittle strength 

• Crazing strength 

• En t an glement molecular weight 

• Distance between entanglements 

• Plateau modulus 

T * >300K ““y ° f the “ P ra P«^ may be calculated as a ftactioa of 

The group parameters required to predict these properties are: 

• Molecular weight 

• Length per repeat unit 

• van der Waals volume 

• Cohesive energy 

• Molar glass transition temperature 
Number of backbone rotatable uni ts 
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QSPR Approaches 


Bicerano’s Synthia Method 

J. Bicerano, Prediction of Polymer Properties (Marcel Dekker, 1993) 


• circumvents reliance on group contributions by relating 
polymer properties to topological information about polymers 
using connectivity indices derived from graph theory 

• useful for any polymer comprised of the following elements: 
C, H, N, O, Si, S, F, Cl, Br 

• applicable to bulk amorphous homopolymers and statistical 
copolymers 


Advantages 

-> independent of large group-contribution library 
-> thermodynamic, mechanical, and transport properties 

Disadvantages 

restricted to nine atom-types 
restricted to amorphous polymers 
-> does no account for special interactions 




Table 11-1. Properties Calculable By the Synthla Module 


Thermophysical Properties 

Electrical, Optical and Magnetic 
ProDerties 

Glass transition temperature 

Refractive index 

Temperature of half decomposition 

Molar refraction 

Coefficient of volumetric thermal expan- 

Dielectric constant 

sion 


Molar volume 

Volume resistivity 

Density 

Diamagnetic susceptibility 

Molar heat capacity at constant pressure 


Cohesive energy 

Chain Stiffness and Entanglement Proper- 

Solubility parameter 

ties 

Surface tension 

Steric hindrance parameter 

Thermal conductivity 

Characteristic ratio 


Molar stiffness function 

Mechanical Properties 

Entanglement molecular weight 

Bulk modulus 

Critical molecular weight 

Shear modulus \ 

Entanglement length 

Young's modulus 


Poisson's ratio 

Transport Properties 

Shear yield stress 

Activation energy for viscous flow 

Brittle fracture stress 

Permeability of gases 
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DESCRIPTION OF METHODS 
Multivariate PLS Approach 
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Multivariate PLS Approach 

Attempts to eliminate any dependence on group contributions or 

connectivity indices 

• relate bulk polymer properties to specific molecular 
descriptors associated with the repeat unit 

• focused on two polymer properties: T g and Tensile Modulus 

• selected seven easily-obtainable molecular descriptors: (1) 
molecular weight; (2) length; (3) number of atoms; (4) number 
of non-H atoms; (5) vdW volume; (6) rotational df, and (7) 
number of backbone atoms 

• develops linear regression equation to map bulk polymer 
property onto calculable descriptors thereby allowing 
prediction of polymer property solely from knowledge of 
descriptors 

property value = a 1 (PCi) + a2 (PC2) + a 3 (PC 3 ) + ... 

• applicable to any series of polymers for which experimental 
data is available 


Advantages 

totally independent of any database library 
generally applicable to any property and any polymer 
-> from Loadings, evaluates contribution by each descriptor to 
property of interest 

Disadvantages 

-> assumes linear relationship between property and 
descriptors 

-4 need expt’l data on related polymers to develop regression 
equation 
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Determining the number of 
significant components 

• cross-validation 
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DESCRIPTION OF METHODS 
Artificial Neural Networks 
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Artificial Neural Networks (ANN) 

paradigm based on architecture of biological network of 
neurons in human brain 


• ANNs learn by repetition (like a child) to generate a desired 
response (output) based on knowledge of specific stimuli 
(inputs) 

• learns to map bulk polymer property (output) onto calculable 
molecular descriptors (inputs), thereby allowing prediction of 
polymer properties solely from knowledge of these descriptors 

• applicable to any series of polymers for which experimental 
data is available 

• focused on same two polymer properties as outputs: T g and 
Tensile Modulus 

• selected identical seven molecular descriptors as inputs: (1) 
molecular weight; (2) length; (3) number of atoms; (4) number 
of non-H atoms; (5) vdW volume; (6) rotational df, and (7) 
number of backbone atoms 


Advantages 

totally independent of any fragment database 
generally applicable to any property and to any polymer 
inherently nonlinear and tolerant of noisy data 

Disadvantages 

requires expt’l data on related polymers to train the network 
prior to making predictions 
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Inputs #1 #2 #3 #4 #5 #6 #7 



T g or E value 


input 

units 


hied-- 


~ \ t p* *• 


output 

neuron 


List of Inputs (per repeat unit) 

#1 molecular weight 

#2 end-to-end length 

#3 number of atoms 

#4 number of non-H atoms 

#5 vdW volume 

#6 rotational degrees of freedom 

#7 number of backbone atoms 
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SUMMARY OF RESULTS 
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Figure 1. Repeat units of polymers included in this QSPR study 
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Table I 

Abbreviated List of Calculated Properties of Polyimides and PEEK 
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Table II. Expanded List of Calculated Properties of Polyimldes and PEEK 
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Table II (continued 



0 Wed on Synthla method. 
















Table III. List of Experimental Properties of Polyimldes and PEEK 
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QSPR-calculated Thermophysical Properties for Various Polyimides 
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TPI, iTPI 



PEEK 


PEEK 

PETI«PCPI*2 
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EXPERIMENTAL POLYMER PROPERTIES 
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CORRELATE WITH EXPERIMENTAL RESULTS 
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Tensile Yield Strength, MPa 


QSPR-calculated Tensile Yield Strength vs. Temperature 
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Poisson's ratio 








ANN and PLS Calculated T g vs. Experimental Values 


Expt’l Calculated/Predicted Tg, K 

T g , K ANN Residual PLS Residual 


CPI 

493 

506 

-13 

473 

20 

TPI 

523 

546 

-23 

520 

3 

Kapton 

658 

634 

24 

664 

-6 

PMR-15 

611 

587 

24 

583 

28 

CPI-2 

490 

513 

-23 

513 

-23 

iTPI 

532 

556 

-24 

550 

-18 

PEEK 

416 

440 

-24 

420 

-4 

PETI-2 


562 


540 


PETI-4 


620 


617 


SI 


519 


633 
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ANN and PLS Calculated Tensile Modulus (E) 


vs. Experimental Values 



EXDt’l 
E, GPa 

ANN 

Calculated/Predicted E. GPa 
Residual PLS 

Residual 

CPI 

4.27 

4.15 

0.12 

4.39 

-0.12 

TPI 

3.7 

3.83 

-0.13 

3.62 

0.08 

Kapton 

2.96 

3.09 

-0.13 

3.32 

-0.36 

PMR-15 

3.9 

3.77 

0.13 

3.65 

0.25 

CPI-2 

3.7 

3.7 

0 

3.53 

0.17 

iTPI 

3.75 

3.62 

0.13 

3.77 

0.02 

PEEK 


3.77 


3.68 


PETI-2 


3.4 


3.49 


PETI-4 


3.97 


4.26 


SI 


4.13 


4.38 
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CURRENT WORK 


• Atomistic Cerius 2 calcuiations of polymer properties: tensile 
modulus, bulk modulus, speed of sound, Poisson’s ratio, 
density 

• EGAMS of neat LaRC-SI resin (with Mike Grayson): 
experimental and theoretical investigation of thermal 
degradation 

• Dynamic mechanical analysis (DMA) of LaRC-SI resin (with 
Mike Grayson): measure T g , characterize sub-T g transitions 

• Extend multivariate regression and neural network analysis 
relating bulk polymer properties to molecular descriptors 
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SUMMARY AND CONCLUSIONS 


QSPR 

• ALL THREE QSPR METHODS AGREE WITH EACH OTHER REASONABLY 
WELL IN PREDICTING MECHANICAL AND THERMOPHYSICAL POLYMER 
PROPERTIES 

• VALUES OF T g AND E PREDICTED BY ALL THREE QSPR METHODS AGREE 
REASONABLY WELL (WITHIN 15%) WITH THE CORRESPONDING EXPT’L 
VALUES 

• COMPARING THE POLYIMIDES, KAPTON RANKS HIGH IN T g AND T m BUT 

LOW IN a T ; Tg IS SLIGHTLY higher for para OVER META catenation 
(488k FOR CPI2 vs. 486k FOR CPI) 

• HIGH Tg AND E ARE INVERSELY RELATED TO a T : KAPTON VS. PEEK 

• LARC-SI EXHIBITS HIGH VALUES OF TENSILE MODULUS AND TENSLE YIELD 
STRENGTH (AND LOW VALUES OF POISSON’S RATIO) ACROSS THE T = 
100-600 K TEMPERATURE RANGE 


PLS REGRESSION AND NEURAL NETS (ANN) 

• FOR Tg, THE RESIDUAL STANDARD DEVIATION (RSD) WAS 1 8K FOR PLS, 
20K FOR ANN, AND 37K FOR QSPR 

• FOR E, THE RSD WAS 0.20 GPA FOR PLS, 0.1 2 GPA FOR ANN, AND 0.91 
GPA FOR QSPR-DOW 

• THE MOLECULAR DESCRIPTOR WITH THE HIGHEST PLS LOADING WAS THE 
ROTATIONAL DEGREES OF FREEDOM FOR BOTH T g AND E 
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RESULTS OF ATOMISTIC 
CERIUS 2 CALCULATIONS 
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Calculated Tensile Modulus : 354 GPa 



Miller planes are highlighted in blue 





University of Akron 



Table C: University of Akron - C Cerius 2 Data 



Enerqv 

density 

Tensile Modulus (GPa) 

Structure 

kcal/mol 

S/cc 

Z axis 

a 

145.0534 

1.1860 

221.3720 

b 

144.7304 

1.1840 

217.1443 

c 

144.972 2 

1.1830 

222.4243 

d 

144.6897 

1.1840 

219.5773 

e 

144.6852 

1.1830 

215.7280 

/ 

144.6677 

1.1840 

219.6206 

g 

145.0247 

1.1860 

219.7376 

h 

145.0469 

1.1860 

200.7676 

/ 

144.7431 

1.1840 

217.9183 

j 

145.0385 

1.1860 

233.4200 

Average 



218.195 
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Chart 1: Predicted Tensile Modulus for Control and Experimental Polymers 



els - PBO Kovlar 49 Akron - A Akron - B Akron - C Akron - D BPDA 

Polymers 





Multiaxial and 
Time DependentResponse 
of 

Rigid Polymers 


W.G. Knauss 


California Institute of Technology 
Pasadena, California 


Supported in part by the Materials Branch at NASA, Langley 
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The following are the captions for the figures. 


1. Tide 

2 . Productions of steel and polymer during last 35 years 

3. Long range objectives 

4. Constitutive Interest 

5. Range of Issues 

6. Yield-like behavior 

7 . Experimental Approaches to Physical Characterization 

8 . A free volume based constitutive model 

9 . Comparison of material responses between the free volume model prediction and 
experimental observation for uniaxial loading and unloading 

10. Comparison of material responses between the free volume model prediction and 
experimental observations for relaxation at different small strain levels 

1 1 . Simulation of yield-like behavior from the free volume model 

12. Free volume model for large uniaxial strains (polymide) 

1 3. Creep compliance at different stress levels measured by Govaert 

14. Master creep compliance curve shifted by Govaert The master curve is not smooth 

15. Comparison of the master creep compliance curve shifted by Govaert and that 
shifted by H. Lu and W.G. Knauss 

16. Zoom of a portion of the master curve "shifted" by Govaert (The curves do not 
really fit a shift) 

17. The shift factors measured by Govaert 

18. The shift factors measured by H. Lu and W.G. Knauss 

19. Displacement measurement resolution of the digital image correlation 

20. Shear relaxation modulus measured at 0.2% of shear strain 

2 1 . Shear relaxation master curve 

22. Young's relaxation modulus measured at 0.2% of axial strain 

23. Young's relaxation master curve 

24. Shift factors for the Young's and shear modulus master curves 

25. Superposed tension or compression increases the shear creep rate under pure 
torsion well below glass transition (results at 22°C) 

26. Superposed tension or compression increases the shear creep rate under pure 
torsion well below glass transition (results at 50°C) 

27. Shear stress also controls the shear creep rate - shear creep compliances at different 
stress levels 

28. Superposed tension or compression increases the shear creep rate under pure 
torsion well below glass transition (results at 80°C) 

29. Under the same torque, a superposed tension increases the torsional creep rate 
while a superposed compression decreases it in the vicinity of the glass transition 
(results at 100°C) 

30. Under the same torque, a superposed tension increases the torsional creep rate 
while a superposed compression decreases it in the vicinity of the glass transition 
(results at 1 10°C) 

3 1 . Under the same shear stress, a superposed tension has a much higher creep rate 
than a superposed compression, indicating free volume plays a major role in the 
glassy state 

32. Some thoughts on an "activation" rate model incorporating the coupling of 
volumetric deformation and shear deformation 
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MILLION TONS 



1960 1980 2000 

YEAR 

Annual Tonnage of Steel and Polymer produced 
during the last 35 years 
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o 

o 


Establish general polymer 
constitutive behavior and 

Analytically model this behavior 


Constitutive behavior must incorporate 

Linear and nonlinear viscoelasticity at 
small and moderate strains 
(yield-like response) 

Nonlinear viscoelasticity at 
very large strains 

under arbitrary load histories, including 
ageing and "recovery" 

Environmental factors: 

Temperature dependence, including the 
range around the glass transition 

Solvent influence, including water 
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Constitutive Interest 


Physical properties descriptions involve the 
whole spectrum of size scales and approaches 
to analytical formulation 


1) Model material behavior for engineering 
analyses: 

Application oriented and global view 

2) Structure - property relationships: 
Molecular science orientation and 
microscopic view 
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Physics-based 
Observation and Issues 


Yield-like behavior is physically different from yielding in 

metals. In thermoplastic solids it is not necessarily 
associated with crazing 

"Yield-like behavior" occurs at small strains (£ U niaxiai<5%) 

"Yield-deformations" are usually recoverable under 
thermal control 

Yield-like behavior does not involve large scale molecular 
re-orientation and organization unless it is 
associated with very localized shear bands 

Molecular orientation occurs in uniaxial deformation at 
extension ratios exceeding 100% 

Very small changes in volume (>0.1%) have a 

disproportionate effect on the time-temperature 
rade-off 

To what degree is yield-like behavior the result of local 

material instability augmented or controlled by small 

volume changes? 

To what degree is yield-like response controlled or 
influenced by organized molecular rearrangement 
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Range of Issues 


1) Time and rate dependence 

2) Nonlinear response 

3) Thermorheology 

4) Glass transition behavior 

5) Physical ageing 
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Experimental Approaches to 
Physical Characterization 


Different stress state histories 


There is a lot of data on pressure superposed on uniaxial 
tension in the literature 

There is no data on biaxial stress states shear-plus-tension 

or shear-plus-compression 


There is virtually no data on time-dependent bulk behavior 
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Constitutive Model (Emri and Knauss, 1981) 


Free volume 

/ = fo + Aa(t) * dT -f 

Time-shift factor, <f> 


lo g<j> = 


b 

2.303 




Internal (reduced) time 


d-m = 


dt 




du 

W(u) ~J{u)} 


Thermorheologically simple material 

Sij=2 IL K( ~ (i) ^dP' dt ' 


— oo 


&kk = 3 
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Log E [bar] 
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M. R. Monaghan 
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Figure 4.6 Three Term Model Tensile Response 
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Figure 4.7 Model Results for Variation in Straining Rates on Specimen of 65mm 

Gauge Length. 
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Displacement Resolution of the Method 


A typical field of view is 25 mm (1 inch) 

Pixel size: 50 microns/pixel side 
The resolution limit for displacements is 3 microns 
This translates into resolution of 
normal strains of 0.05% 
shear strains of 0.05% 
circumferential strains of 0.08% 
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Shear Relaxation Modulus Referred to 105°C 
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Coupling Between 
Volumetric and Shear Response 
as Based on an "Activation" Rate Model 
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Static Properties: RT to (T - 15°C) 
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Micro-Model 
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Micro-Model 


Off-Axis, Test vs Predicted 
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Fracture/Strength RT to (T -15°C) 

• Edge Notched tensile tests 
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Free Volume Evolution During Aging 
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Temperature 




246 




247 




K3B Resin Creep 8-14 
Log Compliance vs Log Time 

Resin Specimen 72-3 Test Temperature 21 5C Stress 387.6 Psi 
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Figure 
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Figure: Comparison between resin and composite long term data at 215°C. 
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Need to expand beyond use of properties from macro 
level materials testing 
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